Geometric features of the mixing of passive scalars at high Schmidt numbers 
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The mixing of passive scalars of decreasing diffusivity, advected in each case by the same three- 
dimensional Navier-Stokes turbulence, is studied. The mixing becomes more isotropic with decreas- 
ing diffusivity. The local flow in the vicinity of steepest negative and positive scalar gradients are 
in general different, and its behavior is studied for various values of the scalar diffusivity. Mixing 
approaches monofractal properties with diminishing diffusivity. We consider these results in the 
context of possible singularities of scalar dissipation. 
PACS: 47.27.Eq, 47.53.+n, 47.27.Gs 



Turbulent mixing of dyes or macromolecular sub- 
stances can be understood by analyzing how these scalar 
substances are advected and diffused in a fluid medium. 
For small concentrations, the admixture, or scalar, is pas- 
sive, or has no dynamic feedback on the flow — which is 
therefore determined independent of the scalar. Often, 
the scalar diffusivity is small compared with the viscos- 
ity of the fluid, so that their ratio, the so-called Schmidt 
number, Sc = v/k, is large. Our interest here lies in the 
advection and diffusion of passive scalars for large values 
of Sc. In particular, we wish to understand the recent 
result from numerical simulations that turbulent mixing 
becomes more isotropic when Sc is increased [1,2]. This 
approach to isotropic state is evident from the fact that 
deviations from it, as quantified by odd moments of the 
scalar derivative along a mean scalar gradient, d^9, de- 
crease as Sc increases from 1 to 64 for a fixed Reynolds 
number of the flow [1]. 

We illustrate the issue on hand in Fig. 1, momentar- 
ily leaving aside the details of how the data were ob- 
tained. The figure shows the probability density func- 
tion (PDF) of d\\9 for four different values of Sc keeping 
the Reynolds number fixed. It is clear from the figure 
that the distribution is quite asymmetric for Sc — 1 
but approaches symmetry as Sc increases. This can be 
seen more quantitatively in the insets. The left inset 
shows the quantity 7r(z) = [p(z) — p(— z)]/[p{z) + p{— z)], 
where z = d\\9 / {(d^) 2 ) 1 / 2 is the normalized scalar gra- 
dient component: smaller values of this quantity indicate 
smaller asymmetries. The right inset shows that a global 
measure of asymmetry A, defined by 



A = 



z 3 p(z)dz/ 



3 p(z)dz . 



(1) 



also decreases with increasing Sc. It is easily shown [3] 
that the asymmetry in the PDF of the scalar gradient is 
a measure of the anisotropy of the small-scale scalar (or 
mixing), so its reduction with Sc is to be regarded as a 
direct indication that mixing becomes isotropic. 

The reason for the asymmetry of the PDF for Sc = 
O(l) has long been known to be the presence of ramp-cliff 
structures in the signature of the scalar field [4] whenever 
there is a mean gradient in the scalar alone or in both 



scalar and velocity. This conclusion has been reinforced 
by recent numerical studies [5-7]. For the present orien- 
tation of mean gradients, the exponential tail on the right 
of Fig. 1 arises directly from cliffs in the scalar signatures, 
the cliffs being related to the fronts of the scalar gradient 
distribution in three-dimensional space. 

Analytical approaches to understanding the features 
shown in Fig. 1 are difficult. Significant progress has been 
made for the so-called Kraichnan model of passive scalars 
[8] in which the advecting velocity field is assumed to be 
Gaussian and to vary infinitely rapidly in time; for a re- 
view see [9] . For two-dimensional and isotropically forced 
high-Sc turbulence it was possible to calculate the PDF 
of the scalar dissipation rate and thus predict the scaling 
in the tails of the PDF of scalar gradient fields [10]. This 
Lagrangian approach was based on a sufficient separation 
of pumping and diffusion time scales which is established 
at large Sc. However, even for the synthetic Kraichnan 
velocity field, a systematic study of anisotropy with Sc 
has not been made (though some beginnings have oc- 
curred with respect to the Sc-dependence in Refs. [11]). 
Another approach [12], which follows Batchelor's original 
model [13] of quasistatic straining motion, provides up- 
per bounds for scalar derivative moments as functions of 
Sc. However, the problem still contains numerous open 
questions. For instance, it is unclear as to what changes 
occur in and around the ramp-cliff structures as Sc — > oo . 
This Letter quantifies some of these changes via numeri- 
cal simulation of turbulent mixing for a range of Schmidt 
numbers. We will particularly consider level sets of the 
steep gradients and relate them to the local structure of 
the advecting flow. The multifractal approach is known 
to sample efficiently the "singular" structures at small 
scales [14-16]. Therefore, we will study the scalar dissi- 
pation field in high-Sc cases and discuss the formation of 
potential singularities. 

Our numerical results were obtained in a homoge- 
neously sheared flow [17] with (u x ) — Cy, where C is 
a constant, at a Taylor-microscale Reynolds number R\ 
of 87 in which the scalar field of constant mean gradi- 
ent was allowed to evolve according to the advection- 
diffusion equation. The Schmidt number was varied from 
1 to 64. The pseudospectral simulations were done with 
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resolutions of 512 x 257 x 512 grid points for a box of 
size 27r : 7r : 27r. We ensure that the scalar fluctuations 
are properly resolved by requiring that k max rjB > 1.3 
with k max = V2N max /3, where N max = 512; the Batch- 
elor scale r\B = Tj/y/Sc and the Kolmogorov scale 77 = 
(^/e) 1 / 4 . The scalar gradient G = e y /ir for all runs. 

To shed some light on Fig. 1, we show in Fig. 2 three 
slices of the scalar field at the same moment of evolution, 
each slice corresponding to a different value of Sc. It is 
clear that the large structure, and the front associated 
with it, do not change with Sc but internal striations 
of finer scale accompany larger Sc. These striations in- 
crease the relative population of negative gradients. This 
is what was observed in Fig. 1: the right side, correspond- 
ing to positive gradients associated with fronts, remains 
effectively unchanged, whereas the left side, associated 
with negative gradients, shows higher probability as Sc 
increases, thus rendering the PDFs increasingly symmet- 
ric. 

The origin of the additional striations for larger Sc can 
be understood by writing down the evolution equation for 
the derivative field, did, where 9 is the scalar fluctuation. 
This equation is 

[d t + u • V - kV 2 ] dfi + d l u ■ V0 = -d,u-G, (2) 

where index i is kept fixed. Compared to the advection- 
diffusion equation for 9, the additional stretching term 
due to the scalar gradient appears. This term (the last 
on the left hand side) is the cause for the amplification 
of the scalar gradient due to local straining motion. 

The authors of Refs. [5,6] related the formation of 
fronts to a straining motion for which the direction of 
compression points into one of the scalar gradients. Sev- 
eral questions arise: How robust is this result for increas- 
ing Scl What type of flow causes the steepest negative 
scalar gradients, and their increase with increasing Sc, 
while the positive gradients remain nearly unaffected? 
To answer these questions, we have performed an eigen- 
value analysis of the local velocity gradients in the vicin- 
ity of the largest scalar gradients. The starting point is 
the characteristic cubic polynomial of the velocity gradi- 
ent matrix 

\ 3 + P \ + q = 0, (3) 

with 

= ' d(u x ,u v ) d(u x ,u z ) d(u y ,u z ) ~ 
[ d(x,y) d(x,z) d(y,z) \ 

and q = —d(u x ,u y ,u z )/d(x,y,z). This equation results 
from det(dui/dxj — A5jj) = 0. The solutions can be di- 
vided into two classes, (i) For A = 27q 2 + Ap 3 > 0, one 
gets a conjugate complex pair of eigenvalues correspond- 
ing to an inward or outward swirl, and one real eigenvalue 
corresponding to the compression or expansion in the 
third principal direction, (ii) For A < 0, all three eigen- 
values are real corresponding to pure strain. For every 
Sc, an interval in the far tails of the scalar gradient PDF 



was used. Figure 3 shows, as functions of Sc, the fractions 
of eigenvalues that were found to be purely real or com- 
plex conjugate; panel (a) corresponds to negative tails of 
the PDF and (b) to positive tails. For positive gradients, 
i.e. for the cliffs, local straining becomes somewhat more 
dominant with growing Sc, while, for the negative tails, 
the swirling and straining motion contribute in almost 
equal parts for large Sc. This is consistent with the view 
that the cliffs are associated with the front stagnation 
point of a moving fluid parcel, where the velocity field is 
predominantly straining, while the negative slopes come 
from the wake region behind such parcels where the fluid 
motion is both vortical and straining. 

Consider now the set F of all scalar concentration gra- 
dients that exceed a certain high threshold value. From 
Fig. 1, we expect the differences between the negative 
and positive gradients to diminish as Sc increases. The 
box counting dimension Do of F, embedded in the three- 
dimensional space, is defined as the scaling exponent of 
N g (F) = N 5- Do where N S (F) is the number of cubes 
of size S that are needed to cover F [18]. The results 
are shown in Fig. 4. While, for Sc = 16, the prefac- 
tor N is significantly different for the two signs of the 
threshold, it is about the same for the two cases for the 
largest Sc, demonstrating that the preference for positive 
gradients — corresponding to the fronts — decreases, and 
that both signs of the gradient appear with increasingly 
equal probability. 

While keeping Sc fixed and varying the threshold we 
observe a decrease of the magnitude of D for increas- 
ing threshold (see lower left inset of Fig. 4). The sets are 
thinned out and their spatial support shrinks to lines, 
i.e. Dq ~ 1. With increasing threshold, the differences 
in Do for level sets with opposite signs increase, which 
indicates that the remnants of asymmetry for even the 
most symmetric case come essentially from the far tails. 

By choosing a range of level sets in the far tails of the 
PDF, \d\\6/G\ > 30tt for Sc = 16, 32 and 64, we ob- 
serve that the box counting dimension decreases weakly 
in magnitude for both signs of gradients (see upper right 
inset of Fig. 4). Based on geometric measure theory [19], 
or just on the constancy of scalar flux [20], it was con- 
jectured that the Hausdorff dimension (which is always 
bounded from above by the box counting dimension) of 
scalar level sets for Sc — > 00 become space-filling. A gra- 
dient level set, which stands for the change from one level 
of the scalar to another, can only have spatial support 
between those of the two levels. If all level sets are space- 
filling, the gradient level sets should show a decreasing 
dimension with increasing Sc. This result does not seem 
to apply for the present instance because we are con- 
sidering here level sets for all gradients exceeding some 
prescribed (though large) threshold. 

The dependence of the box-dimension on the thresh- 
old (lower left inset of Fig. 4) is a clear evidence that the 
scalar derivative is a multifractal. For operational pur- 
poses, it is more convenient to consider the scalar dis- 
sipation field, ee(x,t) — Kj2 3 j = i(dj0) 2 . We define the 
measure 
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(i) 



e#(x + Xj, i)d 3 x 
/v ee(x,i)d 3 x 



(5) 



where Vr^ is the ith cube with length r. The spectrum 
of generalized dimensions [21], D q (q), obeys the following 
scaling relation 



(6) 



In Fig. 5 the results of the multifractal analysis arc shown 
for Sc = 8 and 5c = 64. The scaling range is indicated by 
the dashed lines; though it is small, it exceeds the Batchc- 
lor range and allows the exponents to be determined with 
reasonable confidence. Note that for the largest Sc of 64, 
we have r\jr\B — 8, i.e. less than a decade. The inset shows 
that the D ? -spectrum gets flatter for increasing Sc. We 
have added data points from the dye mixing experiments 
[22] corresponding to a Schmidt number of about 1900. 

The fact that the spectrum of generalized dimensions 
gets flatter with increasing Sc means that in the limit 
of very large Sc the scalar dissipation rate becomes a 
monofractal, i.e. D q — Dq for all q. The physical picture 
is as follows: large "singular" spikes of the scalar dissi- 
pation rate fill out the whole space more and more, so 
that the "quiet" regions in-between become less preva- 
lent. The degree of spatial and temporal intermittency 
decreases, which is just the property to which the spec- 
trum is sensitive. We can quantify this result also by fit- 
ting a bimodal multifractal cascade model to the data 
[23]. For Sc = 8 we obtain p = 0.86, for Sc = 64 a value 
of p = 0.81, decreasing to 0.58 for the experimental data 
at Sc — 1900. The model would yield a monofractal dis- 
sipation field for equidistributed energy flux, i.e. p = 0.5. 
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0121007. 



[1] P. K. Yeung, S. Xu, and K. R. Sreenivasan, Phys. Fluids 

14, 4178 (2002). 
[2] G. Brethouwer, J. C. R. Hunt, and F. T. M. Nieuwstadt, 

J. Fluid Mech. 474, 193 (2003). 
[3] K.R. Sreenivasan, Proc. Roy. Soc. Lond. 434, 165 (1991). 
[4] C. H. Gibson, C. A. Friehe, and S. O. McConnell, 

Phys. Fluids 20, S156 (1977); K. R. Sreenivasan, 

R. A. Antonia, and D. Britz, J. Fluid Mech. 94, 745 

(1979); P. Mestayer, J. Fluid Mech. 125, 475 (1982); 

R. Budwig, S. Tavoularis, and S. Corrsin, J. Fluid Mech. 

153, 441 (1985). 
[5] M. Holzer and E. D. Siggia, Phys. Fluids 6, 1820 (1994); 

ibid 7, 1519 (1995). 
[6] A. Pumir, Phys. Fluids 6, 2118 (1994). 



[7] 

[8] 
[9] 

[10 
[11 

[12 

[13 
[14 

[15 

[16 

[17 
[18 

[19 

[20 
[21 

[22; 

[23 



A. Celani, A. Lanotte, A. Mazzino, and M. Vergassola, 

Phys. Fluids 13, 1768 (2001). 

R. H. Kraichnan, Phys. Fluids 11, 945 (1968). 

G. Falkovich, K. Gaw§dzki, and M. Vergassola, Rev. 

Mod. Phys. 73, 913 (2001). 

M. Chertkov, G. Falkovich, and I. Kolokolov, Phys. Rev. 
Lett. 80, 2121 (1998). 

E. Balkovsky and A. Fouxon, Phys. Rev. E 60, 
4164 (1999); W. E and E. vanden-Eijnden, Physica D 
152-153, 636 (2001); M. Chertkov and V. Lebedev, 
Phys. Rev. Lett. 90, 034501 (2003). 

J. Schumacher, K. R. Sreenivasan, and P. K. Yeung, 

J. Fluid Mech. 479, 221 (2003). 

G. K. Batchelor, J. Fluid Mech. 5, 113 (1959). 

T. C. Halsey, M. H. Jensen, L. P. Kadanoff, I. Procaccia, 

and B. I. Shraiman, Phys. Rev. A 33, 1141 (1986). 

K. R. Sreenivasan and C. Meneveau, Phys. Rev. A 38, 

6287 (1988). 

A. Brandenburg, I. Procaccia, D. Segel, and A. Vincent, 

Phys. Rev. A 46, 4819 (1992). 

J. Schumacher, J. Fluid Mech. 441, 109 (2001). 

K. J. Falconer, The Geometry of Fractal Sets, Cambridge 

University Press, Cambridge, 1986. 

P. Constantin and I. Procaccia, Nonlinearity 7, 1045 
(1994). 

G. L. Eyink, Phys. Rev. E 54, 1497 (1996). 

P. Grassberger, Phys. Lett. A 97, 227 (1983); H. G. E. 

Hentschel and I. Procaccia, Physica D 8, 435 (1983). 

K. R. Sreenivasan and R. R. Prasad, Physica D 38, 322 

(1989). 

C. Meneveau and K. R. Sreenivasan, Phys. Rev. Lett. 59, 
1424 (1987). 




„ 10 

z=3 || e/<(3 || 6) 2 ) 1/2 



20 



3 



FIG. 1. PDFs of the normalized longitudinal scalar deriva- 
tive z for Sc = 1, 4, 16 and 64. The scalar has a mean gradient 
G and is embedded in a Navier-Stokes velocity field, which is 
itself sustained by a constant mean gradient, C . The PDFs 
are shown in the main panel, and line styles in the legend. Left 
inset: the ratio ir(z). Right inset: the global asymmetry factor 
A as a function of Sc. The decay of the derivative skewness, 
(z 3 ), with increasing Sc is comparable to that in [1]. 
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FIG. 2. Slices of the total scalar field (mean plus fluctu- 
ation) for Sc — 1 (left), 8 (middle), and 64 (right), all of 
which are advected by the same flow. Only a fraction of the 
simulation domain is shown. 
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FIG. 4. Box counting for scalar gradient sets F for Sc — 
16 (solid lines) and 64 (dashed lines). The sets are defined 
by |9||0/G| > 3G"7r. Dotted lines indicate the range of scales 
for the fit. Exceeding nominal expectations, the scaling re- 
gion is more extensive than the viscous-convective range (or 
the Batchelor regime). Triangles: positive gradients; asterisks: 
negative gradients. The differences between the number of 
covering boxes for positive and negative gradients are larger 
for lower Sc than for larger Sc. Upper right inset: Box count- 
ing dimensions fitted to the data of the main figure as a func- 
tion of Sc. Lower left inset: Box counting dimension of a scalar 
gradient level set as a function of the threshold value, ci , for 
Sc — 64; here z was taken within an interval c\ ± 0.025. The 
symbols in insets are as in the main figure. 
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FIG. 3. Fractions of eigenvalue triplets, as functions of Sc, 
in the vicinity of: (a) the steepest negative gradients, (b) 
steepest positive gradients. 
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FIG. 5. Multifractal analysis of the scalar dissipation field 
for Sc = 8 and 64. Because all grid points are included, we 
have Do = 3, as indicated by the solid line. The dashed lines 
indicate the scaling range used to determine the D q . The inset 
shows the spectrum of generalized dimensions, D q (q). Sym- 
bols for different Sc are given in the legend. 
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